Quantum Signatures of The Classical Disconnection Border 
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A quantum Heisenberg model with anisotropic coupling and all-to-all interaction has been ana- 
lyzed using the Bose-Einstein statistics. In Ref.pJ the existence of a classical energy disconnection 
border (EDB) in the same kind of models has been demonstrated. We address here the problem to 
find quantum signatures of the EDB. An independent definition of a quantum disconnection bor- 
der, motivated by considerations strictly valid in the quantum regime is given. We also discuss the 
dynamical relevance of the quantum border with respect to quantum magnetic reversal. Contrary 
to the classical case the magnetization can flip even below the EDB through Macroscopic Quantum 
Tunneling. We evaluate the time scale for magnetic reversal from statistical and spectral properties, 
for a small number of particles and in the semiclassical limit. 
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The existence of an energy disconnection border 
(EDB) (previously called non-ergodicity threshold), in 
classical Heisenberg models with anisotropic coupling 
and infinite range of interaction has been recently 
found 1]. Below the EDB the energy surface is discon- 
nected into two components with opposite sign of the 
total magnetization. The dynamical consequences of the 
EDB in classical Heisenberg models with infinite range 
coupling has also been investigated 0). In particular, it 
has been shown that below the EDB the magnetization 
cannot change sign in time, while above it, in a fully 
chaotic regime, a time scale for the magnetic reversal 
(time needed for the total magnetization to flip) can be 
determined. Magnetic reversal times show an exponen- 
tial growth with the number of spins and, as in standard 
phase transitions, a power law divergence at the EDB 
itself. 

The existence of this border is not limited to the infi- 
nite range coupling case and can be in general related to 
the anisotropy of the coupling when it induces an easy- 
axis of magnetization (defined by the direction of the 
magnetization in the minimal energy configuration of the 
system). The relation between the EDB and the range 
of the interaction has also been studied: we have taken 
into account a Heisenberg models with an interaction po- 
tential among the spins which decays as R~ a , where R 
is the distance among the spins. Defining as r, the ratio 
of the disconnected portion of the energy range with re- 
spect to the total energy range, it has been proved that 
for a dimensional system r tends to zero in the ther- 
modynamic limit for a > d(short range) while it remains 
finite for a < d (long range) Q. 

The results found in the classical model guided our 
investigations on the quantum side. We are mainly in- 
terested here in the quantum signature of the classical 
EDB, and on its relevance with respect to the quantum 
reversal time of the magnetic moment. 

We consider here an infinite-range interacting system 
since the explicit expression of the EDB has been ob- 



tained in this case only. Despite its unphysical char- 
acter, magnetic systems can be realized, within modern 
experimental techniques |^ , described by Heisenberg-like 
Hamiltonians with an infinite range term, which could in- 
duce the presence of the EDB. Moreover, when the range 
of the interaction is of the same order of the size of the 
system, the all-to-all coupling could be an important first 
order approximation in the understanding of their behav- 
ior [a 0. This could be the case for small systems used in 
present nano-technology which requires to deal with sys- 
tems with a few dozens of particles 7] , or for macroscopic 
systems with long range interactions. 

We first analyze the spectral properties and we estab- 
lish the existence of a quantum disconnection border in 
close correspondence with the classical one. An analyti- 
cal estimate of this quantum threshold is given. We will 
then study the system from a dynamical point of view, 
analyzing the time scale for quantum magnetic reversal 
and comparing the quantum magnetic reversal times with 
the classical ones. 

We consider a system of N particles of spin I, described 
by the following Hamiltonian, 
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where — 1 < r\ < 1 is the anisotropy constant. We de- 
fine M x , y ,z = J2i S^ ,v,z as the components of the total 
magnetization of the system. Due to the anisotropy of 
the coupling the system has an easy-axis of the mag- 
netization along the y— direction. Quantization of the 
Hamiltonian follows the standard rules. According to 
the correspondence principle, the classical limit is re- 
covered as / — > co. As in the classical case we fix the 
modulus of the spins to one. This can be achieved 
with an appropriate rescaling of the Planck constant, 
h -> h/\Si\ = 1/^1(1 + 1). With this choice, in the 
classical limit, I — > oo (h — > 0), the spin modulus re- 
mains equal to 1. Because of the infinite range nature 
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FIG. 1: (Color online) N = 6l = 3r)=l. a) Doublet 
structure of the low energy region of the spectrum. The dif- 
ferent parity of the states constituting the doublets is shown. 
We also indicated the level distance among even eigenstates, 
A. b) the splitting of the doublets, <5, vs the specific energy, 
e = E/N (blu line) and the nearest neighbor level spacing 
A(e) (red line). Also shown as vertical lines: the quantum dis- 
connection border, e q dis (dashed), computed analytically from 
Eq.@, and edis (dashed-dotted). 



of the interaction, the Hamiltonian P is a completely 
symmetric operator with respect to particle exchange. It 
is thus natural to limit ourselves to subspaces of definite 
symmetry. Specifically, we consider the bosonic case (an 
ensemble of integer spins) , so we will limit our analysis in 
the subspace of all possible completely symmetric states, 
with dimension AT = {N + 2Z)!/(iV!(2Z)!). This choice 
reduces considerably the dimension of the Hilbert space, 
allowing to extend our analysis further in the classical 
limit. An important property of the Hamiltonian Q is 
its invariance under a n rotation about the z-axis : the 
Hamiltonian commutes with the operator exp(i7r^ Sf), 
and its eigenstates can be labeled as odd (-) or even (+) 
according to whether they change or do not change sign 
under such rotation. 

The first aim of our analysis on the quantum system is 
to assess the quantum signature of the classical discon- 
nection threshold, edis- Numerical diagonalization of 
gives rise to a quasi-degenerate energy spectrum with 
a energy splitting increasing with the energy. It is a 
standard result ;8], that the infinite time average of any 
quantum operator is zero in presence of a non-degenerate 
discrete spectrum and if its quantum average over the en- 
ergy eigenstates is zero. The operator M y satisfies this 
conditions thus the total magnetization along the easy- 
axis, can change its sign for any energy. Also, the time 
scale at which this happens can be obtained by a de- 
tailed study of the energy difference between close eigen- 
states. Specifically, since the matrix elements of M y be- 
tween energy eigenstates of the same parity is zero, it 
will be important to study the characteristics of the en- 
ergy distance between even and odd eigenstates. The 
possibility for the magnetization to reverse its sign also 
in the energy region where it would be classically forbid- 
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FIG. 2: (Color online) e q dis (dashed lines) and e* (symbols) 
for different I values vs r\. Data are N = 6, I — 2,4. Also 
shown (as full line) edis that showing, at the same time, the 
agreement between e q dis and e* and their common semiclassi- 
cal limit edis- 



den, can be interpreted as a manifestation of Macroscopic 
Quantum Tunneling j^J , (the total magnetization can be 
a macroscopic quantity) . The evaluation of the tunneling 
rates becomes then crucial to obtain the time scale for 
the quantum magnetic reversal. 

The most evident property of the energy spectrum is 
the presence, in the low energy region, of quasi degen- 
erate doublets, see Fig. Each doublet is composed 
by an even and an odd eigenstate. Even if from Fig. 
they seem degenerate, they are actually split by a small 
energy difference S. At high energy the doublets are not 
well defined anymore. In Fig. QJt we have also indicated 
the level spacing between neighbor even states A. Note 
that the doublets are well defined only when 6 <C A. 

The splittings of the doublets, 5, can be computed nu- 
merically for each even (or odd) eigenstate taking the 
energy difference from its closest odd (or even) state; 5 
can thus be considered as a degree of the spectrum degen- 
eracy. In Fig. \T]p we show how S varies with the specific 
energy, e = E/N: a change of slope is clearly visible, from 
exponential to almost constant. Note that the energy at 
which the slope changes is close to the classical edis: see 
Fig-d 3 • It is possible to have a better estimate of the en- 
ergy at which the slope changes, considering the energy 
e* at which the doublets disappear: 5(e*) ~ A(e*), see 
FigCt>. Note that the value e*, distinguishes two regions 
of the spectrum, characterized or not by quasi-degenerate 
eigenvalues pairs. The behavior of e* for different I values 
as a function of the anisotropic coefficient r\ has been re- 
ported in Fig. |21 As one can see e* — > edis as I increases, 
which confirms the soundness of our definition. 

e* can be thought as a quantum correction to 
the classical EDB. An hand- waving argument allows 
to evaluate this quantum correction: in the classical 
case, edis has been obtained computing the minimum 
of -(V^EOS?) 2 when 77 > 0, and of (r)/2)M% - 

(77/2) J2( s i) 2 when V < °m- Thus > the lowest eigen- 
value of the same operators could give an approximate 
quantum border. We will call this threshold the quantum 
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FIG. 3: (Color online) Average (over energy bins) of (]n(5/D)) 
vs the specific energy e, for N = 6, r\ — 1 and different I, as 
indicated in the legend. Here D is the average level spacing. 
Also shown (dashed and dotted style) the linear fit. Inset : 
the linear dependence of the slopes a with I. Linear fitting is 
q= -1.3 + 61. 
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FIG. 4: (Color online) The quantum average reversal time, 
t (circles), as a function of e is are shown for the case N = 
6, 7} = 1 and I = 4, are compared with the classical ones 
(dashed black line) , showing a good agreement above e q dia and 
a deviation near edis- Also shown, as a full line AP ma x/Po 
averaged over close eigenfunctions. 



disconnection border e q dis , and we have: 

4, ~ -|(«) 2 for V > 

4is ~ l(N-l)(M) 2 for v <0. (2) 

The agreement between the numerical values e* and 
our analytical estimate e q dis has been shown in Fig. [21 
(compare symbols with dashed and dotted lines). 

As shown in Fig. |21 the level splittings increase expo- 
nentially with the energy. Also, on increasing the semi- 
classical parameter I, the rate of growth becomes steep- 
est. In Fig. |3 we show the average of ln<5 over suitable 
energy bins, and normalized to the average level spacing 
D (obtained dividing the energy range by the number of 
states). Linear fits have also been indicated as dashed 
lines. In the inset the exponent a, obtained from the 
fitting : (\n(S/D)) = ae + C has been shown to have 
a linear dependence on the semiclassical parameter I, so 
that 5(E) « for e < edis when I — * oo. 

Let us now analyze the time scale for magnetic rever- 
sal in the quantum system, comparing the results with 
the classical onesQ- Since in the classical case the re- 
versal times have been determined at a fixed energy (mi- 
crocanical approach), we adopt here the same procedure 
and compute the reversal time of the quantum average 
magnetization starting from an ensemble of initial states, 
obtained choosing randomly energy eigenstates in a 
narrow energy interval: \ip) = ^2^ +AE Ce\E). The coef- 
ficients Ce have been randomly chosen in modulus and 
phase and such that ~^2e +AE \Ce\ 2 — 1. Since the total 
magnetization along the easy-axis, M y , connects only 



energy eigenstates with different parity, we have: 

E+AE 

{M y (t)) = 2Ke{ C E+ C E _e-^ T (E + \M y \E_)}, 

E + ,E_ = E 

(3) 

where T = %/ (E- —E + ). ^From (M y (t)) we compute the 
time of first passage through zero for each initial state of 
the ensemble. ^From these times we obtain the average 
magnetic reversal time r. Before presenting the results 
of our analysis let us recall that in the quantum case, at 
variance with the classical one, we are legitimate to ask 
what is the time scale for magnetic reversal in the whole 
energy range. Indeed, since the energy spectrum is non- 
degenerate, from Eq.lO, the average magnetization will 
soon or later reverse its sign, even below the EDB. 

In Fig. 0] we consider the energy region above edis- As 
one can see there is good agreement between classical 
and quantum times above e dis - Note that the classical 
times diverge at edis at variance with the quantum ones 
which are systematically smaller, in the region betwen 
edis and e q dis . This is not surprising since the possibility 
of tunneling will enhance the probability for the magne- 
tization to reverse its sign. In the classical case we suc- 
cessfully evaluated the reversal times from the entropic 
barrier, AS, between the most probable value of the mag- 
netization and its zero value, r ~ e~ AS = P ma x/Po Um- 
Following the classical case we also computed P ma x/Po 
from the quantum probability distribution of the mag- 
netization, P(M y ), and compared the dynamical times r 
with the probabilistic ones given by cP max /Po 0], where 
c is a constant. As one can see from Fig. 0]the agree- 
ment between probabilistic and dynamical times is fairly 
good where classical and quantum times agrees, while 
the agreement is less accurate in the crossover region, 
between edis and e?. . Let us now discuss, in details, 
the behavior of the quantum reversal times below e dis . 
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FIG. 5: (Color online). Comparison between dynamical and 
tunneling splittings in the classically forbidden energy region. 
Circles : dynamical quantum reversal times; Full (blu) line: 
tunneling splittings r ~ nh/(2S) for the case: N = 6, I = 4. 
Not surprisingly, the evaluation of r through S does not work 
above t q dia , where the doublet structure disappear. As dashed 
(red) line P ma x/4:Po- 

In the low energy region of the spectrum, due to the 
intrinsic quasi-degeneracy the dynamics can be entirely 
characterized by the energy difference \E+ — E-\ = 5. 
This occurs if the energy bin AE of the initial state 
is sufficiently small so that one single doublet belongs 
to it. The dynamics is thus oscillatory with a period 
given by 2irh/5. Indeed under this condition, the mag- 
netization oscillates coherently between states with op- 
posite sign, a phenomenon known as Macroscopic Quan- 
tum Coherence 10J. This period also represents, within a 
numerical factor, the time for the first passage to zero of 
(M y (t)). One thus can assume : r ~ tt%/{25). The agree- 
ment, over many orders of magnitude, has been shown in 
Fig- El Also below e dis , we checked the proportionality 
of the reversal times with P ma x/ Po- It is surprising that 
Pmax/Po turns out to be proportional to the tunneling 
rates and then, when properly defined, to the reversal 
times, even in the region classically forbidden (below edis) 
where the only mechanism allowing the jumping of the 
barrier is through Macroscopic Quantum Tunneling, see 
Fig. [3 This suggests that the mechanism producing this 
proportionality can also have a non classical origin. One 
should also note that the constant of proportionality is 
different below edis and above e%, g . This explains the 
poor agreement between the statistical and dynamical 



times in the crossover region (edis < e < e q diB ). 

^From the results presented here we can address the 
problem of the dynamical signature of the classical En- 
ergy Disconnection Border. In the semiclassical limit the 
crossover region becomes very narrow, {e q dis —> £dis)i thus 
we can expect a crossover from of the reversal time from 
a power law dependence on the energy, like in the classi- 
cal case0, to and exponential dependence on the energy. 
Moreover Fig. [3] suggests how the classical limit is recov- 
ered: indeed 6 rj as I — > oo, for energies below edis, 
which is consistent with the fact that the magnetization 
cannot reverse its sign below e d i S in the classical system. 

In conclusion we have found a quantum signature of 
the classical EDB in the spectral properties of the sys- 
tem leading to the definition of a quantum disconnection 
threshold, e dis with the correct classical limit. Below e q dis 
the spectrum is characterized by the presence of quasi de- 
generate doublets, whose energy difference 5(e) depends 
exponentially on e. The quantum reversal times of the to- 
tal magnetization have been studied and compared with 
the classical ones above edis- We have also shown that 
the total magnetization can flip in the energy region clas- 
sically forbidden. Quite surprisingly, quantum reversal 
times (and thus the tunneling rates) are still proportional 
to Pmax/Po even below e dls . 

The existence of the classical EDB allows to ad- 
dress an energy region where to look for Macroscopic 
Quantum Phenomena, which have recently raised much 
interest [llj. Indeed the fact that the total magnetiza- 
tion can reverse its sign even below the EDB can be 
seen as a manifestation of Macroscopic Quantum Tunnel- 
ing, a well known phenomenon in micromagnetism, also 
found experimentally 0, 0]. Nevertheless Macroscopic 
Quantum Tunneling of magnetization arises in literature 
111 llOfl . from phenomenological Single-Spin Hamiltoni- 
ans, where the single spin describes the total magnetic 
moment of the system, and no reference to the range of 
the interaction has been explicitly pointed out. On the 
other hand we presented here a multiparticle system in 
which this phenomenon clearly arises, in connection with 
the existence of the disconnection border and with the 
long range nature of the interaction. 

G.L.Celardo acknowledges financial support from 
LANL and Universita Cattolica under the program For- 
eign specialization studies. 
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